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C^ I evaluated by a Monte Carlo renormalization group method involving a transformation 

"i i in time. The results agree very well with a finite-size scaling analysis performed on 

Q I the same data. The value of z = 2.13 ± 0.01 is obtained, which is consistent with 



most recent estimates. 
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I. INTRODUCTION 

The dynamics of the two-dimensional spin-flip Ising model remains a source of interest 
and controversy. Near the critical point, not only does the correlation length ^ diverge fol- 
lowing C, ~ \T — Tc\~'^ , where Tc is the critical temperature and u a critical exponent, the 
correlation time r also diverges due to critical slowing down. This divergence can be char- 
acterized by the dynamic critical exponent z, where r ~ ^^, or equivalently, r ~ |T — T^^'^, 
where A = uz. The dynamic scaling hypothesis asserts that, in the long time limit, all times 
scale with this diverging time scale, so that a dynamical universality class is characterized 
by, amongst other things, this critical exponent z. Since the particular value of z is common 
to all members of a universality class, it is an important and fundamental quantity. Here 
we use a novel Monte Carlo renormalization group method to estimate z = 2.13 ± 0.01 for 



the two-dimensional spin-flip Ising model, which is in the universality class of model A P9 
where a nonconserved scalar order parameter is the only dynamical mode. 

Many au- 

thors ia|,|,|,|J3J|,|,0,y,|ll,0,0,|l,|Tl,|l3,|ll,0,ig,i3,il,il,il,il,il,i3 have attempted 



to evaluate z by a variety of techniques, giving values that are not always consistent with 
each other, even when considering the given errors. There are many possible explanations 
of these discrepancies, and we shall consider some of them. First, however, we shall discuss 
some lower bounds. For model A, mean fleld theory predicts z = 2 — 1], where t] is the critical 
exponent describing the power-law decay of the correlation function. Using a generalized 
Langevin-equation approach, Schneider |^ showed this was a lower bound on z. Of course, 
the mean-fleld result is correct at and above the upper critical dimension du = 4, below which 
the e = 4 - rf expansion |2| gives z = 2 + 0.01345(4 - d)^ - 0.002268(4 - df + C(4 - d^. 
The lower critical dimension of model A is (i^ = 1, and Bausch et al. |2^ have calculated the 
d= l-|-e expansion of the kinetic drumhead model to be: z = 2+{d—l) — ^{d—l)'^+0{d~l)^. 
More recently, one of us [Q argued that the dynamic critical exponent should be larger than 



the reciprocal of the exponent for domain growth, yielding z > 2. 



Of the different techniques used for the determination of z, no method seems to have 
proved better than others. Moreover, a given method does not always yield consistent re- 
sults. For example, high-temperature expansions of different orders apparently converge to 
different values whether the expansion is computed for spin-spin time correlations |2^J2^ or 



magnetization correlations pB|,p7|. On the other hand, Monte Carlo techniques have sys- 



tematic errors which can be quite difficult to evaluate. For example, it is generally believed 
|17| J3T|JT5|J3^ that the dynamic universality class of the spin-flip Ising model is insensitive to 



the algorithm used as long as the updating algorithm is local [l33| , |3^ . Some results obtained 
in one dimension [^,^ support this hypothesis. However, studies in three dimensions 



using Creutz's deterministic microcanonical dynamics |^ yield a value of z that "agrees 
with the Monte Carlo measurements as they agree among themselves" but is nevertheless 



slightly higher than most recent estimates (e.g. |37| , |38[| and references therein). This dis- 
crepancy could be due to the methods of estimation of z rather than to the algorithm, but 
to our knowledge, there exists no systematic analysis testing the limits of validity of the 
dynamic universality class hypothesis. This is a potentially important issue, since nonlocal 
acceleration algorithms are being developed [§^,0 which have exceedingly small values of 
z. 

Besides the algorithm, another source of systematic error is from critical slowing down 
itself. While such effects are well known and are typically incorporated |]39| , ^ to consistently 



estimate 2, the existence of long time correlations can potentially induce subtle couplings of 
correlations in pseudo-random number generators. 



It has also been noted l2^:^l|] that the linear relaxation critical exponent of the magne- 
tization M can differ from the nonlinear one. The nonlinear relaxation time r is defined 
as r = /Q°°(M(t) — M(oo))/(M(0) — M{oo))dt where the denominator is for normalization. 
If the integrand is written as a sum of exponentials, J2i Wie~*''^^ , then r will simply be the 
weighted sum of the partial relaxation times. The limit when M(0) — > M^q , i.e. , as the 
system is left to relax from a value near its equilibrium value, defines the linear relaxation 
time. It characterizes a relaxation process that does not include the scaling of the order 



parameter itself. Racz et al. p3| , p]| used scaling arguments to suggest that the nonlinear 
critical exponent A)^ of an observable A can be related to the linear critical exponent A^ 
by the following relation Aj^ = A^ — jSa where I3a describes the scaling of the quantity A 
with respect to temperature. Mori et al. [^ in d = 2, and Chakrabarti et al. |^2| in three 



dimensions, reported observations of such a relation. 

Current techniques used in Monte Carlo simulations sometimes involve nonlinear re- 



sponse in terms of the definition given above. Some authors []T2|,|ll],|16|,|T8[ have used the 
relaxation of the order parameter in systems prepared at zero temperature when put in 
contact with a heat bath at some temperature near T^. On the other hand, other techniques 
||l|j^j5|,|^,[f|,p|,p|, p!0| , p!2| , pr3| involve the measurement of some time correlations in critically equili- 
brated systems. Presumably there could be differences involving transient dynamics in these 
methods. However, there is no real distinction made in the literature concerning possible 
discrepancies between these various methods. It is not clear yet how important the dis- 
tinction between linear and non-linear relaxation times is, and further investigations should 
clarify this situation. 

Finite-size effects can be important in simulations of critical systems, e.g., the effect of 
using the infinite system critical temperature can induce systematic errors in small systems. 
Although these effects can be exploited by finite-size scaling, very small systems may not 
be in the range where scaling applies. Finite systems also introduce the concept of ergodic 



time |[T7| , |43[| , i.e., the mean lifetime of the system in one of its broken symmetry states. 
While infinite systems at T^ have a vanishing order parameter, finite systems of size L 
have nonzero values ±|Mj^| between which the system has spontaneous transitions. Those 
transitions introduce a large effect in time correlations as discussed below. 

Finally, it should be noted that experiments have been done on systems thought to be 



in the universality class of model A. The few results [0,0 of which we are aware yielded 
values of z well below 2. However the experiments are difficult, and the microscopic processes 
involved in the systems are various, so that it may be that the experiments do not probe 
the problem of interest herein. 



II. METHOD 



The Monte Carlo renormalization group (MCRG) was introduced by Ma |i6], and devel- 



oped by others, especially Swendsen [^]. For critical dynamics, the method was extended 
by Tobochnik, Sarkar, and Cordery |^, and others [^,12,18,^. MCRG allows one to use 
the self-similarity in critically equilibrated systems by analyzing the effect of a controlled 
change of length scales on correlation functions. When length scales are changed by a factor 
of b, by some suitable blocking of b"^ spins to one renormalized spin, the correlation length 
changed hj C, -^ C,/b. This implies that time scales are changed by a factor of r — > r/b^, 
from which the dynamic critical exponent can be estimated. Unfortunately the method 
is self-consistent in that there is no proof the system approaches a fixed point under the 
renormalization group. Thus it is essential that checks are made that, after several levels of 
KG, scaling is consistently observed. Here, we generalize the usual procedure of blocking in 
space, by blocking in time t. 

Consider the Hamiltonian of an Ising-like system 

n = J2K^s^ (1) 

a 

where the K^s are the coupling constants including temperature, the a index runs over all 
i = 1,2,- ■ ■ ,N spins for nearest-neighbor, next-nearest neighbor (and so on) interactions, 
and the SaS are generalized spins made of specific products of spins a^ = ±1 on each site. 
For example, the Ising model has Ki = —J/kBT, and Ka>i = where J is the coupling 
constant, ks is Boltzmann's constant and T is temperature and Si = a^aj such that the sum 
is restricted so that i and j are nearest neighbors. Here we will consider the two-dimensional 
Ising model on a square lattice, and apply a renormalizat ion-group transformation repeatedly 
to this evolving system. As mentioned above, a typical numerical MCRG transformation 
46| is to "block" by a length rescaling factor of b. A block of b''' spins is transformed into a 



renormalized spin by majority rule of the spins in the block, with a random outcome on ties. 
The resulting renormalized Hamiltonian is assumed to be expressible in terms of another 



short-range Ising-like system, with more Sa terms contributing. Numerically, the approxi- 
mation results from the fact that the number of spins remaining after m blockings is N/b"^'^, 
thus coupling constants for interactions are truncated, if they involve larger length scales 
than the entire system. Thus, for equilibrium properties, a suitably chosen transformation 
operator will change the Hamiltonian Ti^K) by moving the parameter vector K to some 
other point in 1^-space. 

For critical dynamics the situation is analogous. Say one begins with the Ising model 
again, with the dynamics of single-spin flips, where each evolving configuration depends only 
on the previous one, i.e., a Markov process. Under the blocking transformation mentioned 
above, it is expected that the original master equation which is "short-range" in time is 
changed to a non-Markovian equation which has memory over some small time scales. 

Our original contribution here is to introduce blocking in time as well as space. We 
simulate a process on a regular time scale measured in terms of Monte Carlo steps (MCS). 
In addition to applying a standard real-space MCRG transformation technique, we perform 
a blocking of spins, by majority rule, in consecutive discrete time steps. In principle, the 
advantage of blocking in time can be two-fold. Firstly, we expect that it will smooth out high 
frequency fluctuations allowing one to reach the asymptotic limit of the RG transformation 
more rapidly than with blocking in space alone. In the same way as blocking in space iterates 
away irrelevant short-length-scale behavior, we expect that blocking in time will further 
eliminate short-time irrelevant memory effects. Secondly, one can suitably adjust the time 
blocking factor ht to balance the effects resulting from blocking in space, by choosing, for 
example, ht = b^°, where Zo is some reference exponent. We choose Zo = 2 for reasons that 
will become clear below. 

Before discussing the method in more detail, we introduce the quantities that will be 
measured in time. Critical dynamics involves a scaling relation in which all time scales, in 
the long time limit, are measured in units of the diverging correlation time. Thus, one needs 
only to choose a convenient measure of time correlations. We shall use the time displaced 
correlation function for the magnetization M, as defined by 
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(. {MjQMjU + t)) - {M{Q){M{U + t)) 

'^''^ ' {Mito) - {M{Qfy/^M{t, + t) - {M{to + t)fY/^' ^ ' 

where M is Yh c"i as usual. Other correlation measures are discussed below. All averages 
were computed from selected numerical discrete time steps 5t (usually a few MCS). In 
principle, the time-correlation function should be fitted to a series of exponentials 

VM{t) = Y.^^e-"^' (3) 

i 

where r, > Tj if z < j. However, tests made to fit the data to two exponentials showed that 
Ti ^ T2. Indeed, all our data were well fit to a single exponential with time constant tm- 
Data were extracted from one long simulation from which values of v^m(^) were computed 
over a time range of t = to t = 5rj\//. Further averaging was also made by running 32 
systems in parallel. 

To calculate z^ we use a matching procedure |§,PJ^. In principle, after the irrelevant 
variables have been iterated away, the probability distribution function will remain invariant 
under further renormalization-group transformations. It is expected that, after a finite 
number of iterations, contributions from the irrelevant variables will be negligible. Then, 
any quantity determined after m blockings of an N spin system should be identical to those 
determined after m + 1 blockings of a system of Nh'^ spins. In our new method, time scales 
after m blockings of the small lattice are also explicitly rescaled by a factor {bt)"^, while the 
larger lattice has times explicitly rescaled by a factor (64)'""*"^. Unless bt = b^, quantities 
measured on the two lattices will still be at different times t and t'. Hence, close to the 
fixed point, we expect a matching condition to hold: !f{N,m,t) = ip{Nb'^,m + l,t') for a 
correlation function (f. From this, the time rescaling factor t' /t can be calculated, through 
the measurement of a suitable correlation function, ipM in our case. The difference between 
the estimate of z and the correct dynamical critical exponent can then be obtained, since 

^ = fe^-^°, (4) 

where b^° = bf. 



Our RG transformation was done in the following way. During a simulation, every four 
configurations, each separated by one MCS, were "blocked" in space and time: One block 
spin was made by majority rule of the 16 spins coming from 4 consecutive configurations of 4 
neighboring spins, ties were broken at random. This corresponds to a space blocking factor 
6 = 2, and a time blocking factor of bt = 2^° = 4. Such a choice would give asymptotically 
trivial rescaling of length and time if 2; = 2, and we expect it to make our study sensitive to 
the difference {z — 2), which is small. 

Instead of doing point to point matching, i.e., matching each discrete time step, the 
quality of our data is such that we have matched the entire function ipM{t), since it could be 
well fitted to one exponential. More explicitly, we have tm ~ ^^ and a renormalized system 
for which r|^ ~ ^'^ where ^' = ^/b. Therefore, without blocking in time, the critical exponent 
is obtained from tm/t'j^j = V ■ Now, if time is rescaled in such a way that r'l^ = r'jy.f/bf, then 
t'Ii/tm = V /bt = b^~^°. The discrepancy {z — 2) can then be obtained from 

\n{TMiL, m)) - ln{TM{bL, m + 1)) 

'-^ = hTb • (^) 

Simulations were done on two-dimensional nearest-neighbor square lattice spin-flip Ising 
systems with periodic boundary conditions. We used a single-flip multi-system algorithm 
i.e, one running different systems in parallel instead of the more common multi-flipping 
one-system algorithms. The only correlation between the parallel systems is the sharing of 
the updating sequence history but we expect this to be negligible. The dynamic algorithm 
was of Metropolis type, i.e., one using a flipping probability based on min(l,exp~^'^/^^^). 
Systems were initialized for 20-50 r^/ at the critical temperature of the infinite system [^ . 
All measurements were made at the same temperature. 

The results could also be interpreted using a finite-size scaling analysis. According to 
this approach, the correlation length of the systems should be of the order of the system 
size so that 

tm-^L^ (6) 



Therefore, the correlation time for different system sizes directly yields the critical dynamic 
exponent. Consistent results from both methods will be presented in the next section. 

III. RESULTS 

From each of the simulations detailed in table | a correlation time was extracted, as 
listed in table f|. Figures |^ and ^ show a typical decay of ipM(t), demonstrating that the 
time correlations for finite systems can be well described by an exponential. All the fits 
were done using a least-square fit algorithm. The dynamic critical exponent could then be 
obtained by comparing the values of tm for different systems. Values of z thus obtained are 



listed in table |T|. The first striking fact is that a point to point finite-size scaling yields a 
value that compares very well with the one obtained from the dynamic MCRG we propose. 
This suggests that the discrepancies in the values of z as obtained from systems of different 
sizes are not systematic errors in the evaluation methods. 

Although our study involves more accumulated averages than any previous work, it is 
still difficult for us to uniquely extract errors. We have found that the "instantaneous" values 
of Tm contains large fluctuations. Figure § shows that even when binned and averaged over 
~ 75tm, the value of tm averaged for 32 systems still contains a fair amount of fluctuations. 



Therefore, as is well known |p9| , ^ , very long simulations are required to get a representative 
value of time correlations. Furthermore, the regimes of fitting of v^a/ (^) have been chosen as 
well as possible, but it should be noted that a logarithmic scale applied to small numbers 
can introduce large fluctuations. Nevertheless, we found that our data was well represented 
by simple exponential decay over the long times we considered (t < 5tm)- 

If one computes the same time correlations with the absolute value of the magnetization, 
then one finds much smaller correlation times, thus showing that the ergodic time is for the 
most part responsible for the large value of tm- Since the ergodic time is due to a finite-size 
effect, it is not surprising to find such a good agreement between our results and finite-size 
scaling. Also, comparison with other kinds of time correlations M, e.g. 
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c(t) = Y.^M<^o + t) 



or 



ij 

where i and j are nearest neighbors, showed that (pMit) is by far less noisy. See figure 0. 

Evaluation of the data with standard finite-size scaling gives a result of 2; = 2.14 ± 0.01, 
consistent with recent works which used this technique [^,|Ty,|I^] on systems of comparable 
size. However, pointwise scaling shows that small systems tend to show a larger value of z 
and that this tendency seems to disappear for L larger than ~ 32. This result is somewhat 
surprising since a previous MCRG study observed a small z ^ 2 from matching systems of 
size L = 8 and 16, and obtained a larger value when matching larger systems. On the other 



hand, a numerical study of very large systems |T^ gave a low value oi z = 2.076 ± 0.005. 
Therefore, there is still some reason for concern on z^s dependence on system size. 

Some authors [^ have argued that the smallness of the dynamic critical region could 
be held responsible for the difficulties encountered in evaluating z. However, the scaling 



relations obtained near the critical point |]T^ do not seem to show such a sharp region, and 
those effects are thought to be strong finite-size effects. Moreover, our MCRG results are 
quite consistent from the first level of RG, as the level of iteration of RG is increased. This 
self-consistently shows that the system is in the critical regime. 

Finally, by considering that small systems have a systematic error and do not seem to be 
in the asymptotic scaling regime, we obtain an estimate of 2; = 2.13 ± 0.01 from the MCRG 
method we propose. This value is consistent with most of the recent work P,|18|,^,||,|10|,[12|,|13 



as well as our finite-size scaling study herein. 

IV. CONCLUSION 

We showed that time can be used as a renormalizable variable in a MCRG method. 
Combined with standard real-space MCRG techniques, we extracted a value for the dynamic 
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critical exponent that was consistent with values extracted from the same data by finite-size 
scaling. Our estimated value oi z = 2.13 ±0.01 is also consistent with most recent estimates. 
The motivation for introducing this method was the possibility of obtaining a more accu- 
rate value of z. Unfortunately, we found that the method was not superior to conventional 
finite-size scaling or real-space MCRG techniques for this problem, although it proved to 
be at least equivalent. Further studies would be required to see if this new MCRG method 
would be of use for other related dynamic problems, such as spinodal decomposition. 
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FIGURES 

FIG. 1. 

Different time-displaced correlation functions for a system of 64 x 64 sites. Note the large 
difference between <fE{t) and ^pM{t) showing that the energy relaxes much more rapidly than 
the order parameter. Also note the noise common to C{t) and E{t). We used <^m(^) in our 
estimations. 

FIG. 2. 

Critical dynamics MCRG on a 64 x 64 system. Curves are for m = 0,1,2,3,4, from 
top to bottom. Averaged over 32 independent systems observed for 11 468 800 mcs. The 
equilibrating time of 204 800 mcs and ipM calculated every 16 mcs. 

FIG. 3. 

Critical dynamics MCRG on a 32 x 32 system. Curves are for m = 0, 1, 2, 3, from 
top to bottom. Averaged over 32 independent systems each observed for 8 192 000 mcs. 
Equilibrating time of 81 920 mcs and ifM calculated every 16 mcs. 

FIG. 4. 

The evolution of the value of the cumulated average of the time-displaced correlation 
functions for a system of 32 x 32 sites. Even if we accumulated a large amount of data, we 
see that the cumulative average still contains a fair amount of fluctuations. 
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TABLES 
TABLE L Simulation details for the various systems. The range is the amount of Monte Carlo 

steps per spin (mcs) for which the function (puit) has been extracted. It can be thought of as 

an observation window over one simulation running in time. The total mcs for one system can be 

obtained by multiplying columns 2 and 4. All the systems were first equilibrated for 10 times the 

value of column 2, which in turns is of the order of 2-5 r. 



L 


range (mcs) 


5t (mcs) 


average # 


w-max 


8 


1024 


4 


32000 


2 


12 


2048 


8 


32000 


2 


16 


8192 


16 


32000 


3 


16 


8192 


16 


32000 


3 


16 


4096 


16 


32000 


2 


20 


8192 


16 


32000 


2 


24 


8192 


16 


26400 


3 


24 


16384 


16 


32000 


3 


32 


16384 


16 


32000 


4 


32 


16384 


16 


32000 


4 


48 


32768 


16 


4832 


4 


64 


65536 


64 


18016 


4 


96 


98304 


64 


2336 


4 



TABLE II. The values of tm as estimated from a least square fit of the time-time correlation 
function (fMit) to a simple exponential. The values are in mcs and the simulation characteristics 
can be read from the respective entry in table 1. The errors indicated in parenthesis are those 
obtained from the fit. The last line represents the value of z obtained from finite-size scaling 
analysis applied to each column. 



L 


MCRG iteration number 







1 


2 


3 


4 


8 


261.41(6) 


64.31(4) 


16.34(4) 






12 


631.8(2) 


157.9(1) 


39.5(1) 






16 


1193.9(5) 


298.4(2) 


74.45(4) 


18.6(1) 




16 


1179.3(3) 


294.6(2) 


73.5(1) 


18.5(1) 




16 


1189.1(5) 


296.9(3) 


74.3(2) 






20 


1930.7(3) 


482.9(1) 


120.5(2) 






24 


2853.9(3) 


716.0(4) 


179.0(3) 


44.8(2) 




24 


2837.9(3) 


713.3(6) 


178.3(3) 


44.5(1) 




32 


5366(2) 


1341.9(6) 


335.4(3) 


84.0(2) 


20.7(2) 


32 


5343.4(5) 


1335.6(3) 


333.8(2) 


83.4(1) 


20.8(1) 


48 


12528(2) 


3130.4(5) 


782.8(3) 


195.6(1) 


49.4(2) 


64 


23316(4) 


5814(2) 


1453(1) 


363.3(8) 


92(1) 


96 


54620(20) 


13660(10) 


3415(5) 


852(3) 


214(2) 


z 


2.153 


2.145 


2.142 


2.125 


2.128 
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TABLE III. Values of z as estimated from point to point finite-size scaling (FSS) and by 
matching the correlation times according to relation (^). 



L 


MCRG iteration number 




FSS 


1 


2 


3 


4 


8-16 


2.183 


2.182 


2.182 


2.182 




12-24 


2.171 


2.178 


2.178 


2.179 




16-32 


2.172 


2.173 


2.174 


2.176 


2.167 


24-48 


2.139 


2.138 


2.132 


2.131 


2.144 


32-64 


2.124 


2.120 


2.118 


2.119 


2.137 


48-96 


2.124 


2.125 


2.126 


2.122 


2.130 
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Figure 1 
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Figure 2 
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Figure 3 
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Figure 4 
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